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ABSTRACT 

We have derived the bivariate luminosity function for the far ultraviolet (1530A) and far infrared 
(60 /xm). We used matched GALEX and IRAS data, and redshifts from NED and PSC-z. We have 
derived a total star formation luminosity function 4>(L to t), with L to t — Lpjjv + Lpm. Using these, 
we determined the cosmic "star formation rate" function and density for the local universe. The total 
SFR function 4>(L to t) is fit very well by a log-normal distribution over five decades of luminosity. 
We find that the bivariate luminosity function 4>(Lpuy, Lpj R ) shows a bimodal behavior, with Lpm 
tracking Lpuy for Ltot < 10 10 L©, and Lpijy saturating at ~ 10 10 L Q , while Ltot ~ Lpi R for 
higher luminosities. We also calculate the SFR density and compare it to other measurements. 
Subject headings: Ultraviolet: galaxies Infrared: galaxies galaxies: fundamental parameters galaxies: 
luminosity function, mass function galaxies: evolution 
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1. INTRODUCTION 

The evolution of the cosmic star formation rate (SFR) 
density represents a fundamental const raint on the 
growth of stellar mass in galaxie s over time ijMadau et al.1 
ll996UFall. Chariot, fc Peil9 961. The distribution of star 
formation rates, or the "SFR Function" , in galaxies is po- 
tentially also a fundamental constraint on cosmological 
models and on the physics of star formation in galaxies. 

While a number of SFR metrics have been used in the 
past, perhaps the most direct measurement of SFR is the 
bolometric luminosity of massive stars, usually obtained 
from the sum of far ultraviolet and far infrared luminosi- 
ties. With the launch of GALEX, a large, homogeneous, 
magnitude limited sample of UV measurements can be 
combined with the IRAS FIR sample to generate a true 
bolometric luminosity function in the local universe. As 
discussed in this volume by Bu at et al.l ((2004), FUV and 
FIR selected samples have quite distinct far ultraviolet 
(FUV; 1350-1750A) to far infrared (FIR; 60^m) lumi- 
nosity ratios. However, if the samples are large, homo- 
geneous, flux limited, and deep enough in both bands, 
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volume dependence can be removed and the fundamen- 
tal bivariate distribution derived for either sample. To 
provide th e most information, the samples can also be 
combined l|Avni and BahcaHlll98flD . 

Our goal in this paper is to generate a bolometric lu- 
minosity function and luminosity density for the bands 
which sample recent star formation. We do this by us- 
ing FUV-selected, FIR-selected, and combined samples 
to estimate the bivariate luminosity function (BVLF) 
in L FUV and L FIR , (f>(L FUV ,L FIR ), using the V max 
method. We then bin this BVLF into a single, to- 
tal luminosity function (TLF) (f>(L to t), where L to t = 
Lfuv + Lfib- We provide some simple parametric fits 
for the TLF, and estimate the cosmic SFR density. We 
discuss whether FUV and FIR selected samples provide 
consistent measurements of these functions. We conclude 
with a brief discussion of the implications. Our cosmol- 
ogy is Q m = 0.3, Q A = 0.7, H = 70 km/s/Mpc. 

2. SAMPLES 

We used two samples to generate the TLF: a far-UV- 
selected sample (FUVS) and a far-IR-selected sample 
(FIRS). The FUVS consists of a primary FUV-selected 
sample and a IRAS FIR matched co-sample. The FIRS 
consists of a primary FIR-selected sample and a GALEX 
FUV-selected matched co-sample. We measured aper- 
ture fluxes for all co-sample matches using optical catalog 
ellipses. A small fraction in each co-sample are formally 
non-detections-the effect of including these (negligible) 
is discussed in §3. 

The FUVS was generated using GALEX All-sky Imag- 
ing S urvey (AIS) and Medium Imaging Survey (MIS) 
data ijMartin et a II l200l iMorrissev et all l200|) . The 
sample consists of objects in GALEX Internal Release 
0.2 (consisting of 649 AIS and 94 MIS pointings) with 
FUV<17 that have a catalog entry in NED. The FUV 
magnitude limit was selected to insure that, if the 
galaxy was not detected in the FIR, an IRAS SCANPI 
ijHeiou. Khan. Malek. fc Boehmerlll988tlBeichman et alJ 
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1988) upper limit (~0.1 Jy) would be meaningful. NED 
was used to determine the galaxy redshift and size. An 
elliptical aperture based on the galaxy size (usually from 
RC3) was used to determine the FUV magnitude, since 
the patchy nature of galaxies in the FUV occasionally 
leads to object shredding by the GALEX pipeline. We 
verified that the requirement for a NED entry did not 
compromise completeness or the FUV-selected nature of 
the sample. We did this using the GALEX/SDSS-DR2 
overlap. Out of 32 FUV<17 objects that SDSS classi- 
fied as galaxies, 31 objects were in NED (the other was 
a close degenerate star). 

The IRAS 60 and 100 micron data for the FUVS 
has been compiled from the following sources - listed 
in th e order of preference. 1. Bright Galaxy C at- 
alog ijSoifer. Boehmer. Neueeb auer 1 fe Sanders! I1989|) 2. 
IRAS Large Opti cal Catalog iRice et alJll988D 3. Faint 
Source Catalog iMoshir fc et all 119901) . 4. SCANPI 
(v2.4). Of the 220 galaxies, 83 had no published 
IRAS fluxes. Detections >3 sigma were extracted from 
SCANPI processing for 37 of the 83. The remaining 46 
are 3 sigma upper limits (all are >0.1 Jy) as measured 
from the SCANPI processing. With regard to SCANPI, 
the median coadded scan flux values are always used. We 
assume all are extended sources. 

The F IRS was generated u sing PSC-z as the primary 
catalog (jSaunders et al.ll2000j) . GALEX FUV data from 
the AIS was available (as of 4/8/2004) for 3938 deg 2 
of the all-sky PSC-z catalog. In the overlap area, 991 
galaxies appear in the PSC-z catalog, with 878 having 
valid redshift and FUV data, for a completeness of 89%. 
This small incompleteness is unlikely to affect the re- 
sults. Two methods were used to correct for shredding 
of la rge galaxies, both using the APM ellipse parame- 
ters llMaddox. Efstathiou. Sutherland, fc Lovedavlll990t 
iSaunders et alJ |1990l) . which axe based on second- 
moment fitting and are scaled to equal the APM detec- 
tion isophotal area (24 mag arcsec -2 for O-plates and 
23 mag arcsec~ 2 for E-plates). In the first method, all 
GALEX catalog objects found within the optical APM 
ellipse were summed to produce a total magnitude. In 
the second, an aperture magnitude was obtained within 
the APM ellipse multiplied by two. The latter also pro- 
vides FUV fluxes for PSC-z objects with no GALEX- 
detected FUV counterpart (112 out of 878 objects). As 
we discuss below, including these non-detections (our de- 
fault) does not affect the results. 

3. LUMINOSITY FUNCTIONS 

We calculate the bivariate and total luminos ity func- 
tions using the 1/V max weighting method JSchmidt 
1968). V m ax is calculated for each object 0Villmcr 
1997D . and for the FUV and the FIR limits of the sam- 
ple. For the FUV limits, we account for the field ex- 
posure time and local extinction (standard GALEX cat- 
alog value). A higher extinction reduces V max [FUV]. 
The adopted V max depends on the treatment of non- 
detections. For both FIRS and FUVS, we created 
samples that excluded the non-detections and included 
them. When excluding non-detections, we use V max = 
min(V m ax{FUV},V max [FIR}), since an object can only 
be detected in both samples if it falls within both vol- 
umes. When including non-detections, all sources are for- 
mally included with flux estimates obtained in the identi- 



cal fashion as true detections. In this case, V max for the 
co-sample is formally infinite. For the FUVS, V max = 
V max [FUV], and for the the FIRS V max = V max [FIR}. 
Finally, both the BVLF and the TLF are obtained by 
summing 1/V m ax into logarithmic luminosity bins, with 
AlogL = 0.5. 

For simplicity in this preliminary study, luminosi- 
ties are defined as follows: Lpuv = ^fuvLfuv,^, 
where Lfuv.u is the monochromatic FUV luminosity; 
Lfir = vqoLqq.v, where L^o^ is the monochromatic 
luminosity at 60 /im. The total luminosity is defined 
as Ltot — Lfuv + Lfir, and converted to star for- 
mation rate usin g SFR[M© yr _1 ]= 3.5 x 10 9 Ltot [L ] 
l)Kennicuttill998|) . More complex relations (e.g., Lfir — 
0. 65^60-^60, i/+0. 42z^ioo£ioo, i/) produce similar results but 
with more dispersion with respect to the simple func- 
tional fits discussed below. We make no k-corrections, as 
rcdshifts are quite low (z < 0.04). 

The TLF has been calculated for the FUVS (168 ob- 
jects including non-detections, 136 excluding them), the 
FIRS (878 objects including non-detections, 766 other- 
wise), and a combined sample. We find that the results 
do not depend on the the inclusion of non-detections, so 
all results we present include them. If all three sam- 
ples include representative galaxies present in the lo- 
cal universe, the derived luminosity functions should be 
consistent within errors from sampling and cosmic vari- 
ance. The combined sample is gener ated following the 
"incoh erent" combination method of lAvni and Bahcalll 
(1980). The FUVS is obtained in a subset of regions cov- 
ered by the FIRS. Thus the FUVS adds information only 
for objects that have f$Q < 0.6, the PSC-z 60 /im limit. 
The combined sample therefore consists of the FIRS sam- 
ple plus the FUVS (/ 60 < 0.6) subsample (113), for a 
total of 991 objects. 

The three TLF are displayed in Figure Q As hoped, 
the three samples give quite consistent results. The 
FUVS slightly exceeds the FIRS at lower luminosities, 
while the FIRS fills out the high luminosity end not rep- 
resented in the fairly bright cutoff FUVS. We restrict our 
analysis here to the combined sample. The error bars are 
generated using a bootstrap method. 

As is apparent in Figure ^ a Schechter function pro- 
vides a very poor fit to the higher luminosity portion of 
the TLF, giving a total x 2 = 11 5 for 9 d.o.f. On the 
other hand, a log-normal function ijSaunders et alJll990(l 

j, /n ,, T </>* -(logL/L*) 2 

4>{L)d\ogL = — 7H= ex P 7T- 2 \d\o%L (1) 

a v 2tt 2eH 

provides a remarkably good fit. In this case, the pa- 
rameters 0* = 0.150 ± 0.035, logi, = 7.43 ± 0.17, and 
a = 0.87 ±0.03 (errors generated by the bootstrap) yield 
a total x 2 — 8 for 9 d.o.f. 

The luminosity distribution L<j){L) is also log-normal, 
peaking at L'„ = L* + a 2 In 10, or L'„ = 9.37, as we show 
in the inset of Figure ^ R can be seen that 50% of the 
SFR density comes from galaxies with log Ltot < 9.4, 
or about 1 M & /yr. 

We have estimated the BVLF using the combined sam- 
ple, and we show a 2D histogram normalized by the TLF 
in Figure [2^. The histogram shows quite dramatically 
why the local FUV LF is well fit by a Sche chter func- 
tion ijWvder et all 120041: llrever et al. Il2004|) : the FUV 
luminosity appears to "saturate" at Lpuv > W 10 Lq, 
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with all increase in the total luminosity coming from FIR 
radiation. This saturation apparently occurs at higher 
redshift, but at a factor of 20 higher Lpuv f° r z =3 
iJAdelbereer and Steidel2000() . For L FUV < 1O 1O L , the 
FIR and FUV luminosities track, but with a slope steeper 
than unity. There appears to be a trough between these 
regions. The trough is statistically significant as it falls in 
the range in which the object number distribution peaks. 

In Figure we show the BVLF rebinned in logarith- 
mic Ltot and Lfir/ Lpuv bins. The total luminosity 
is well correlated with the FIR/UV ratio, and a line is 
shown with the fit Ltot = 9.4 + 1.3(log Lfir/Lfuv) — 
0.15 * (log L fir/ 'Lpuv) 2 ■ The trend of increasing 
LfjrI Lfuv with increa sing Ltot was first noted by 
Wang & Hcckmanl lITool . 

Note also that the FUV projection of the BVLF for the 
FUVS sample is in excellent agreement with lWvder et all 
(2003). 

We use the BVLF to calculate the luminosity density 
using a simple sum: 

p l = J f L l (f)(L FUV ,L FIR )dlogL FUV d\ogL F iR. From 
the inset in Figure ^ it is apparent that extrapolation 
to low or high luminosity using a model fit to calculate 
the luminosity density would not alter t his result greatly. 
The results, using the lKennicuttl itlQOl SFR conversion 
factor, are [L FU v, L F ir, L TO t] = [0.010±0.0014, 0.011± 
0.0005, 0.021 ± O.OOl^MQyr^Mpc- 3 . Hence the lumi- 
nosity density is split roughly 50/50 into primary FUV 
and reprocessed FIR light. 

4. DISCUSSION 

We have made the first attempt at deriving the bivari- 
ate luminosity function for the two bands which trace the 
high mass star formation rate in galaxies. The BVLF 
can be derived from FUV, FIR, or combined samples. 
The resulting functions agree for the samples we stud- 
ied, as hoped. Large, homogeneous, combined samples 
that probe the bulk of the BVLF will provide an excel- 
lent tool for studying the relationship between FUV and 
FIR emission. 

We have used the BVLF to generate a total high 
mass star formation luminosity function and luminos- 
ity density for the local universe. Our value for the 
star formation rate density, 0.021 M^yr^ 1 Mpc~ 3 (< 
z >~ 0.02), is in good agre ement with the estimate of 
LPerez-Gonzalez et alJ ()2003[) using extinction corrected 
Ha of 0.025 M i/r _1 Mpc~ 3 a nd that using extin ction 
corrected FUV from GALEX l|Wvder et al]l2004tK also 
0.025 M Q yr- 1 Mpc- 3 . 

There are a number of striking features of the BVLF 
and TLF. We have pointed out the divergence of Lfir 
and LpuVi behavior similar to that observed when 
comparing FIR and opt ical light ijSaunders et aljfl990t 
iBuat fc Burgarellalll99^ . The BVLF has a bimodal ap- 
pearance, roughly divided at 3 Mgyr^ 1 . Perhaps below 
the threshold SFR, star formation is an equilibrium pro- 
cess and feedback is successful at clearing sightlines for 
FUV emergent flux. With a very high SFR, the process 
may take on a non-equilibrium character, where feedback 
fails to (or has yet to) clear paths for primary radiation. 
The FIR/FUV flux ratio appears to be a rough total 
luminosity proxy. One explanation is that the highest 



SFRs occur in the most massive star forming galaxies 
ijBrinchmann et alJ2003D. The most massiv e galaxies are 
the most metal rich l|Tremonti et alJl2004l . and a high 
metallicty ISM has a high gas to dust ratio. Also, the 
highest SFRs appear to occ ur in galaxies with the highest 
ISM surface mass densities l)Kennicuttll 1989^1 . and higher 
dust column density. 

It has also been known for some time that the 
FIR luminosity function was not a Schechter func- 
tion ((Saunders et al.lll990(l . whereas we no w know from 
GAL EX that FUV luminosity function is ijWvder et alJ 
2004). Taken together with the diversity of star forma- 
tion modalities making up our sample, it is therefore sur- 
prising to learn how well the total luminosity function 
is described by a sin gle log-normal functi on, over five 
decades of luminosity. iNorman et all (|2004f) found a log- 
normal distribution for the X-ray luminosity function, 
and notes that it is the expected distribution for a com- 
plex multiplicative random process. Perhaps a unifying 
physical framework can be found for star formation in 
galaxies ranging from irregular dwarfs to ultra-luminous 
merging galaxies. 

Optical and near IR band luminosity functions that 
trace stellar mass are well fit by Schechte r func- 
tions (Bell. M cintosh. Katz. fc Weinberg! |2003|) . as is 
the UVLF. The fact th at the TLF is so distinct 
from that of the fuel HI llZwaan et alJl2003h and CO 
ijKeres. Yun. fc Young|l2003|) . which are also well fit by 
Schechter functions, implies that the star formation effi- 
ciency is much higher in luminous star forming galaxies, 
a widely accepted result. The fact that the TLF rises 
far above the stellar mass LF argues strongly that the 
timescale for star formation is a falling function of the 
SFR, perhaps due to fuel consumption or feedback. The 
BVLF in Figure 3a is clearly a major clue to the ultimate 
origin of the high luminosity cutoff of the Schechter func- 
tion. 

With this background, it is extremely interesting to 
study how the BVLF evolves over cosmic time. Ly- 
man break galaxies (LBG) show rest FUV well-fit by 
Schechter functions, but with a characteristic UV lumi- 
nosity a factor of 20 larger. While there is no question 
that luminous star forming galaxies were more pervasive 
in the past, a major question in extragalactic astron- 
omy remains the relationship between galaxies selected 
by rest UV vs. rest FIR. Our work suggests that with the 
launch of GALEX and the Spitzer Space Telescope, and 
the flowering of Sub-mm astronomy, a unified approach 
to combining rest FUV and FIR information will bring 
major insights in the next few years. 
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Fig. 1. — The total star formation luminosity function for the local universe. Black points show the FUVS, red points the FIRS, and 
green points the combined sample, for the samples that include non-detections. Error bars are shown for the combined sample. A histogram 
shows the number of galaxies in each one-half decade luminosity bin. The low luminosity bin with one galaxy has been suppressed. The 
curve is the best fit log-normal function with 4>* = 0.15, log L* = 7.43, and a = 0.87. The histogram shows the number of FUVS (blue), and 
FIRS (red) sources in the combined sample in each luminosity bin, with numbers indicating the total in each bin. INSET: The luminosity 
density distribution function Ltot4>{Ltot) for the combined sample, with a line showing the theoretical distribution using the TLF 
log-normal model fit parameters. 
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Fig. 2. — LEFT: a) The bivariatc luminosity function 4>(Lpuy, Lpm), normalized by 4>(Ltot to compress the dynamic range. RIGHT: 
b) The BVLF rebinned as </>' (Ltot)i ^firI Lfuv)> again normalized by 4>(Ltot)- The line shows a quadratic fit discussed in the text. 
In both cases the grayscale density scales linearly with the normalized distribution. 



